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In this article, we analytically study second harmonic (SH) generation 
from thin metallic films with subwavelength, non-centrosymmetry patterns. 
Because the thickness of the film is much smaller than the SH wavelength, 
retardation effects are negligible. The far-field SH intensities are thus domi- 
nated by an effective electric dipole. These analytical observations are further 
justified numerically by studying the effect of polarization of the fundamental 
field on both the SH signal and the electric dipole. It is demonstrated that 
bulk SH polarization density is comparable with its surface counterpart. The 
electric dipole, consequently, originates from the entire volume of the metallic 
membrane, in contrast to the fact that SH generation from metal surface is 
generally dominated by a surface dipole. © 2009 Optical Society of America 
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Optical second-harmonic generation (SHG) from one-dimensional metal surfaces were first 
observed in 1965 [1], and attracted considerable attentions in the following fifty years (see 
Ref. [2, 3] and the cited references). Recently, scientific interests have gradually shifted 
to the quadratic nonlinearities of membrane-like metal films with subwavelength, non- 
centrosymmetry patterns, partially owing to the significant near-field enhancement induced 
by the excitation of localized plasmon resonances [4-6]. On the experimental side, SHG 
was observed from different geometric configurations such as split-ring resonators [7, 8] and 
their complementary counterparts [9], noncentrosymmetric T-shaped nanodimers [10], T- 
shaped [8] and L-shaped nanoparticles [11,12]. 

For ideally infinite metal surfaces, it is well known that the dominant second-harmonic 
(SH) electric dipole source appears only at the interface between centrosymmetric media 
where the inversion symmetry is broken. Higher order multipole sources, such as magnetic 
dipole and electric quadrupole, merely provide a relatively small bulk SH polarization density. 
In contrast to an infinite surface, we will show in this article that there are volume electric 
dipoles which dominate SHG from thin metal films with non-centrosymmetry patterns. The 
underlying physical reasons include: (1) the non-centrosymmetry patterns allow SH electric 
dipoles to appear in the overall structure [13, 14]; and (2) the thickness of the film is far 
smaller than the SH wavelength thus retardation effects are negligible. 

In general, SH radiation of a patterned structure satisfies the following inhomogeneous 
wave equation [15], 

V X V X B{u) - — E = zc^^o/(r)j(r, u), (1) 

where /(r) stresses the structural geometry: it equals 1 for metal and for vacuum. We now 
consider an ideal current sheet whose thickness (along z direction) rfj, <^ A with A being the 
SH wavelength. The SH electric field at direction n in the far zone is found to be [16] 

E = '^/''i^ X p) X n, (2) 



with T] = ^/Jlo/eo and 

p{n) = J J i{rii,z)e-"'^''r-\\+^^^)(mdz. (3) 

Here k = 2'k/X is the wave number in vacuum, and the integration is performed over the 
volume of the sheet. At either positive or negative z direction, we have 

The far-zone field is therefore transverse with a vanishing z component. Without loss of gen- 
erality, we choose the electric field along the x direction (namely, p\\ = Pwe^)- The associated 
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SH intensity is thus proportional to 



g{z)e """dz 



at +2; direction; 
a.t —z direction, 



where the integral function 



9{z) = J j^{rii,z)drii 



(5) 



(6) 



stands for the total current at a specific z plane. 

To identify the contributions of different order multipole sources, we expand Eq.(5) in a 
power series of kdz [15]. At the leading order, the intensities at the ±2 directions are identical 
and proportional to 



g{z)dz 




j^{Y\\,z)dY\\dz 



(7) 



The right-hand side of Eq.(7) is the total current of the sheet or equivalently the effective 
electric dipole. At the first order of kd^, these two signals have a difference of 

g{zi)g*{z2){zi - Z2)dzidz2, (8) 



2ik 

lo Jo 

and the associated relative difference reads as 



Im 



-Ak J^^ g{z)zdz 
it 9{z)dz 



(9) 



In other words, the difference is induced by the appearance of magnetic dipole and electric 
quadrupole. Their influence is negligible under the conditions such as (1) complex g{z) with 
slowly-varying phase, or (2) dz/X -C 1, or (3) g{z) is symmetrical such that g{z) = g{dz — z); 
or (4) g{z) is localized at one point (such as SHG from metal surface). On the other hand, 
g{z) with rapid-varying phase, such as a nonlocal electric dipole, can lead to a remarkable 
difference [13, 17, 18]. However, we will show it is unlikely to happen in an ideally thin film. 

To numerically study second-order nonlinearities of thin metallic films, we employ a clas- 
sical model developed recently [3,9]. This model has been demonstrated to not only provide 
qualitative agreement with experiments, but also reproduce the overall strength of the exper- 
imentally observed SH signals. Before we present the numerical results, it should be stressed 
that, as suggested by simulations (not shown here), these results are quite general and can be 
found from different configurations including (1) perfect T-shaped, C-shaped and L-shaped 
patterns; (2) ideal Ci^-symmetric patterns supported by semi-infinite dielectric substrate; (3) 
asymmetric pattern with perturbed Ci^-symmetry; and (4) particle arrays with size-tolerance 
disorder. 
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The representative film studied here is a freestanding array of gold split-ring resonators 
with a thickness of 25 nm, which closely matches the experimental sample in Ref. [7-9]. 
Schematic drawing of this patterned film is shown in Figure 1, together with its linear 
spectra. In the following, we choose a fundamental-frequency (FF) field which propagates 
along the z direction and has a wavelength of 1212 nm. This wavelength corresponds to the 
valley in the transmission which is induced by the fundamental plasmonic resonance [19]. 

Figure 2 shows the effect of the incident polarization angle Q on both the far-zone SH 
electric field as well as the total SH current. The distinct dependence between the different 
component of the SH electric field can be interpreted in terms of the electric dipole approxi- 
mation [3]. The most striking feature is that almost identical SH electric fields, with a relative 
difference as small as 0.5%, are found along the (±2;) directions. Moreover, the amplitude 
of the electric dipole is proportional to the SH electric fields. Consequently, SH radiation 
from this thin film is dominated by the electric dipole defined in Eq.(7), and higher-order 
multipoles can be safely neglected. It should be mentioned that the same observations are 
also experimentally found from an individual metal tip: the emission of SH radiation at the 
tip can be attributed to a single on-axis oscillating dipole [20,21]. 

To explore the characteristics of the electric dipole, the spatial distribution of jy with a 
x-polarized illumination is calculated (6' = therefore = 0) and plotted in Figure 3, at a 
time corresponding to a maximal jy. We further define the following functions 



with zq = 12.5 nm. Here g{z') stands for the total current at the z = z' plane, and l{x') 
represents the sum of jy along the line (x = x',z = zq). gs{z) is a cumulative function with 
Qsizmax) Corresponding to the electric dipole. Clearly, the g{z) here is almost symmetrical 
in 2, mainly because the film thickness is so thin that both FF and SH fields are nearly 
constant inside the film. In other words, the retardation effects are negligible. As discussed 
previously, it is this symmetrical g{z) which results in the negligible higher-order multipoles. 
On the other hand, each element of the metal volume contributes to the SH radiation, in 
a constructive/destructive way. Assuming the outermost grid layer (a finite-difference time- 
domain method is utilized here [3]) constitutes the metal surface (its 2.5-nm thickness is 
much thicker than a real surface), then the internal bulk polarization density is found to 
be comparable with its surface counterpart. As a consequence, the effective electric dipole 
originates from the entire metal volume. In addition, SH current distribution on the z = 12.5 
nm plane is shown in Fig. (3c). Around the inner corners, the current is found to be strongly 
depolarized and form localized bright hot spots. This phenomenon has been observed from 




(10) 



y 



x'<x 
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metal surfaces with nanoscale roughness, and the origin is attributed to the overlap between 
SH and fundamental modes [22,23]. 

We now discuss the experiment reported in Ref. [11] where considerable multipole in- 
terferences, in addition to the dominant electric dipole, were observed. The experimental 
sample is a 20-nm-thick gold film periodically perforated with an L-shaped pattern, and the 
FF wavelength is 1060 nm. The phase variation of the FF wave inside the metal is small, 
approximately O.OItt, which should result in a g{z) quite close to the one shown in Fig. 
(3a). Consequently, the multipole contribution should be negligible. The difference between 
experiment and theory is believed to be induced by sample imperfections which may result 
in strongly localized out-of-phase currents at both (according to the condition (4) mentioned 
above) top and bottom metal-dielectric interfaces. On the other hand, we want to point 
out that significant multipole sources may appear in large-size metal particles [13, 17] and 
multilayered structures such as bulk photonic metamaterials [24]. 

The authors thank Stephan W. Koch and Colm Dineen for their helpful comments. This 
work is supported by the Air Force Office of Scientific Research (AFOSR), under Grant 
No. FA9550-07-1-0010 and FA9550-04-1-0213. J. V. Moloney acknowledges support from the 
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Fig. 1. (a) Schematic drawing of a square lattice of gold split-ring resonators. The thickness 
of the gold film is 25 nm. (b) Its linear spectra under a x-polarized normal incidence. 
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Fig. 2. The effect of incident polarization angle 9 on the far- zone second harmonic electric 
field along the ±z direction, with ^ = corresponding to polarization along the x direction. 
The fitted function employed is (a) cos^ 9 and (b) | sin 26' |, respectively. The averaged current 
density (proportional to the electric dipole) j is also shown. The amplitude of the incident 
fundamental frequency electric field i?o is 2 x 10^(V/m), same as that used in the experiments 
[7]. The bulk plasma frequency of gold is taken as uOp = 1.367x lO^^s"^, the phenomenological 
collision frequency 7 = 6.478 x lO^^s"^ [3]. 
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Fig. 3. The spatial distribution of second-harmonic current (with arbitrary units) with a 
x-polarized normal incidence. The resolution is limited to 2.5 nm due to computational 
resource limitations. z{x) = corresponds to the starting plane of the gold film along the 
z{x) direction, (c) The current distribution on the plane of Zq = 12.5 nm. 
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